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ABSTRACT 

Stars may be assembled in large growth spurts, however the evidence for this hypothesis 
is circumstantial. Directly studying the accretion at the earliest phases of stellar growth is 
challenging because young stars are deeply embedded in optically thick envelopes, which have 
spectral energy distributions that peak in the far-IR, where observations are difficult. In this 
paper, we consider the feasibility of detecting accretion outbursts from these younger stars by 
investigating the timescales for how the protostellar envelope responds to changes in the emission 
properties of the central source. The envelope heats up in response to an outburst, brightening 
at all wavelengths and with the emission peak moving to shorter wavelengths. The timescale 
for this change depends on the time for dust grains to heat and re-emit photons and the time 
required for the energy to escape the inner, optically-thick portion of the envelope. We find that 
the dust response time is much shorter than the photon propagation time and thus the timescale 
over which the emission varies is set by time delays imposed by geometry. These times are hours 
to days near the peak of the spectral energy distribution and weeks to months in the sub-mm. 
The ideal location to quickly detect continuum variability is therefore in the mid- to far-IR, near 
the peak of the spectral energy distribution, where the change in emission amplitude is largest. 
Searching for variability in sub-mm continuum emission is also feasible, though with a longer 
time separation and a weaker relationship between the amount of detected emission amplitude 
and change in central source luminosity. Such observations would constrain accretion histories 
of protostars and would help to trace the disk/envelope instabilities that lead to stellar growth. 

Subject headings: 



Introduction 



A st ar accretes most of its mass during a ~ 0.5 My r period while it is deeply embedded within a dense 
envelope ([Evans et all 120091 ) . In the seminal IShul (|1977l ) model of spherical protostellar collapse, the young 
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star grows steadily fr om the infalling env elope, with the fixed mass-infall rate determined solely by the 
envelope sound speed. iKenvon et al. found, however, that the luminosities of most protostars fall far 

below those expected from energy release by steady accretion over protostellar lifetimes. They suggested 
that rare a ccretion events ( outbu rsts) could reconcile this discrepancy. Using large photometric samples of 
protostars, iDunham et al.l (|2010l ) confirmed the severity of the "luminosity problem", and found that the 
observed luminosities are well explained when invoking accretion outbursts, which may account for ~ 90% 
of the total growth in stellar mass. The protostellar luminosity fun ction may instead be explained by 
competitive accretion , which increases the duration of weak accretion ( Offner fc McKeel 2011 : Mversl 2011 : 
Krvukova et al.ll2012h . 



The episodic accretion hypothesis has significant support from later stages of pre-main sequence stellar 
evolution. Spectacular outbursts of FUors a nd EXors are usual l y interpreted as episodes w i th accretion 
rates enhanced by factors of 100-1000 (e.g., Herbig 1977 . 20081 : Hartmann fc Kenvon 1996 : Aspin et al 



2010f ). Models of disks suggest that such outb ursts may be triggered by the combina tion of the gravitational 



instability and magneto-rotation al instability (lArmitage et alj|2001t IZhu et al.ll2010 ), although orbital evolu 



tion of binary star/disk s ystems (IReipurth &; Aspinll2004l ). gravitational clumping (|Vorobvov fc Basull2005l) . 



and thermal instabilities (IBell et al.lll995l ) could also explain these violent outbursts. Disk instabilities are 
likely not limited to only those large outbursts of accretion. For classical T Tauri stars, accretion is also 
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2012t ICostigan et al.ll2012j) . Such variability is also expected for younger, more embedded sources. 



Despite many well-known cases of large accretion outbursts seen at later evolutionary stages, they are 
difficult to directly stud y in the youngest protostars. Searches for outbursts are easiest t o perform in the 
optical and near-IR (e.g. McLaughlin! 1946 : Covev et alJboil : Morales-Calderon et alJlioil ). Although some 
FUors and EXors retain remnant envelopes, most of the stellar growth is thought to occur when the central 
protostar is too deeply embedded to detect in the optical or near-IR. Evidence for outbursts in the youngest 
protostars is th erefore mostly indirec t, with the strongest support from bolometric luminosity-temperature 
diagrams (e.g., IDunham et alJ 120101 ). Indeed, measured accretion rates from embedded objects are not 
significantly larger than t hose from CTTSs and are ~ 2 orders of magnitu de too low to explain the growth 
of a star in 10 5 — 10 6 yr (jCaratti O Garatti et alj|2012t ISalvk et al.ll2013l ). Additional support for episodic 
accretion can be found in periodic shocks that occur along the axis of protostellar jets , which may provide 
a record of the past accretion history of stars (e.g. Reiourth 19891 Raga et al. 2002 ). and from chemical 
abundances of Stage objects (|Kim et alj|201lt IVisser fc Berginll2012al ). 



Models of periodic outbursts are able to explain the distribution of bolomet ric temperatures and lumi 



nosities for protostars (jDunham et al.N2010t lOffner fc McKeel 1201 It iMversI 1201 II ). although models of time 
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of periodic outbursts can be tested directly by finding and characterizing these explosive events during the 
initial stages of stellar growth. Moreover, if the accretion process onto the youngest protostars is similar to 
that onto older ones then similar disk instabilities or changes in the star-disk interaction should also lead to 
frequent changes in their accretion luminosities. Measuring the frequency and scale of larg e accretion events 



will help to assess the physics of large-scale instabilities in the disk (e.g. IZhu et al] 120101 ) . These measure- 
ments are also necessary to determine whether accretion bursts are large enough to reduce the timescale for 
stellar contraction, w hich reduces the pre-main sequence lifetime and may introdu ce scatter into pre-main 



sequence stellar ages ( Baraffe et al. 2009 : Hosokawa et al. 2011 : Baraffe et al. 20121 ) 



However, the dense envelope that surrounds the youngest protostar clouds whether any variability in the 
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central source would even be detectable. The physical and chemical envelope models described above assume 
that after an accretion outburst, the envelope reaches a new equilibrium within their model timcsteps, often 
~ lOOyr. The accretion luminosity from the central protostar is reprocessed by dust in the envelope into 
far-IR and sub-mm photons. Our ability to directly detect any change in the accretion luminosity from 
a growing protostar depends on the average travel time for energy to get absorbed by dust grains in the 
envelope and then re-emitted by dust at longer wavelengths, eventually escaping. In this paper, we discuss 
the prospects for detecting such changes by analyzing time-dependence in the radiative transfer of photons 
in the envelope of a deeply embedded protostar. 



2. Observing Envelope Variability 

In order to investigate the manner in which time variability of accretion onto a deeply embedded 
protostar affects the brightness and spectral energy distribution (SED), wc begin with a simple envelope 
model and consider both the steady-state emission profile and the timescales required for achieving steady- 
state. 



2.1. Fiducial Envelope Model 



The IShul (|1977l ) inside-out collapse theory for the evolution of a protostar surrounded by an isothermal 
spherical envelope provides a simple semi-analytic model. The initial density distribution p s in the static 
envelope is power-law in form, with 



Ps 



2nG 



(1) 



where c s (T) is the isothermal sound speed in the envelope, T is the temperature of the gas in the envelope, 
G is the gravitational constant, and r is the radial distance. At time t x from the start of the inside-out 
collapse, material inside of R x = c s t x will be infalling with a shallower density power-law such that over 
most of the infall it has the form 



3 n z 

o=- s B~ 1/2 r- 3/2 
At the same time, the central protostellar source will have accumulated a mass 

,2 



5 

G~ 



given the constant infall rate 
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We further consider the initial envelope mass to be M cnv , which we attain by truncating the outer portion 
of the singular isothermal sphere. Thus the outer radius of the static envelope is 

G 



R e 
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(5) 



For our fiducial model we need to define three parameters, (M env ,t x ,T), or alternatively (M env , M ps , T). 

The luminosity L ps produced by the protostar at the cen t re of the envelope can be estimated from its 
mass M ps and an appeal to protostellar models such as lStahlerl (jl983l) . For accreting low-mass protostars the 
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intcrnal radiation is augmented by the accretion luminosity L acc , which should be significant and is given as 

_ GMpgMps , , 

^acc — 5 , \0) 

-ftps 



where the protostellar radius is also tabulated bv IStahlerl (11983). Although epi dosic accretion may decrease 



the photospheric luminosity relative to the steady state (e.g. lBaraffe et al.ll2009j ). this change from the steady 
state would only reduce the age of the star in our fiducial model. Combined, the steady-state luminosity 
from the central source is expected to be 

Ltot = L acc + L ps . (7) 

The emergent spectrum of the source will depend upon the temperature of the photosphere and the contri- 
bution and shape from the accretion flow. In the spherical envelope case, however, the extremely optically 
thick inner protostellar envelope effectively reprocesses all of this radiation. Thus, the emergent SED should 
be independent of the SED from the central rcgiorQ. 

The opacity within the envelope is provided by dust. We assume that the dust opacity follows that 
for thin icy mantles mixed within a gas density of 10 6 gcm -3 (commonly referred to as OH5 opacities 



Ossenkopf fc Hennina 119941 ) and that the dust mass is one percent of the gas mass. The optical depth 
through the envelope is dominated by the material residing at small distances from the protostar. Thus 
understanding the morphology of the envelope on small scales is important. The envelope is expected to 
have an inner cavity surrounding the protostar due to the presence of an outflow (see discussion in §3.3). As 
well, dust close to the protostar will evaporate when the temperature is high, removing the primary source 
of opacity in the most inner regions of the envelope. 

For the spherical envelope model, the exact location of this inner opacity cavity does not play a significant 
role in determining the emergent spectrum. Most of the energy from the protostar is initially emitted as 
optical photons and then is reprocessed to much longer wavelengths before escaping the system. In this case 
the envelope can be effectively broken into two pieces: (1) an inner zone where the effective temperature 
in the envelope is such that the phot ons produced are trapped, and (2) an outer zone where the photons 



produced are free to escape. Following iHartmannl (|1998l ). the approximate location of the boundary between 
these zones, i? p h, can be determined by requiring that at this radius the temperature be approximately that 
of a black body radiating such that L tot —An R^ h o~ T 4 while, at the same time, the optical depth outward 
through the envelope from this location for photons at the peak of the blackbody distribution is r ~ 1. We 
refer this location as the effective photosphere of the envelope. 

We take as our fiducial model a 1.5 Mq envelope with an isothermal temperature of 10 K, resulting in 
a total size of i? env = 1-8 x 10 4 AU. The embedded protostar is taken to have a mass of 0.25 M Q and a 
luminosity L ps = 1.2 L Q , reflecting an age of about 1.6 x 10 5 yrs after the start of collapse. Steady-state 
accretion should produce £ acc ~ 5 L Q but as noted in the introduction, observations of deeply embedded 
protostars do not in general reveal such high luminosities, implying lower steady-state accretion rates. In 
order to test the importance of variability in the accretion rate, we consider three scenarios: L to t = 1-2, 12.0, 
and 120. L , referring to them as Quiescent (no or little accretion, equivalent to a rate Ox steady-state), 
10 x Burst (enhanced accretion, equivalent to a rate 2.2 x steady-state), and 100 x Burst (very enhanced 



1 We note that while the protostar will heat the inner region of the envelope to temperatures much higher than the isothermal 
value used to compute the initial density distribution, the heated gas is infalling and has little opportunity to react in response to 
this physical change. Thus, the derived isothermal density distribution should hold even though the gas is no longer isothermal. 



- 5 - 



accretion, equivalent to a rate 12 x steady-state), respectively. The values chosen are only meant to be 
illustrative of the underlying processes while the theoretical determinations below show the expected scaling 
for any burst luminosity. For ease of reference, the physical properties related to the fiducial model are 
tabulated in Tabic [TJ 

Wc use DUSTY \lvezic fc Elitzull997l) to compute the equilibrium dust temperature profile throughout 



the envelope (Figure [T]) for each of these cases as well as the emergent SEDs (Figure [5]). No external 
interstellar radiation field is included in these envelope models, so the outer envelope can get unrealistically 
cold. We return to this issue in the Discussion, Three aspects of the plots are immediately obvious. 
First, the higher the central luminosity, the hotter the envelope at all radii, with the temperature increase 
roughly consistent with T oc L 1 / 4 . Second, the higher the central luminosity, the brighter the SED at all 
wavelengths. Finally, the higher the central luminosity, the shorter the peak wavelength o f emission. For an 



excellent diffusion approximation explanation of these consequences see iHartmannl (|1998l ) 



Figure [2] suggests that the most logical place to look for evidence of protostellar variability is near the 
peak of the SED, at mid- to far-IR wavelengths, where the change in the brightness of the embedded protostar 
varies almost linearly with the change in the source luminosity. For significantly shorter wavelengths the 
reprocessing in the extremely optically thick interior masks any ability to observe the protostellar variation 
(although this might change dramatically for non-sphcrical envelopes as we discuss in i} |3.3l bclow). At longer 
(millimetre) wavelengths the spectrum is dominated by the Rayleigh- Jeans tail of the emission from the outer 
envelope, where the bulk of the envelope mass resides. Thus the increase in the emitted radiation at long 
wavelengths varies approximately linearly with the increase in the outer envelope temperature (see Figure 
[T]). However, the observability of variations in the envelope at long wavelengths is worse than it appears in 
Figure ED since the outer parts of real protostellar envelopes are heated by the external interstellar radiation 
field and therefore reach a roughly constant equilibrium temperature, independent of the internal luminosity. 
We return to this issue in the Discussion, §3.11 

The peak of the SED in Figure ED varies from about 100 urn for the Quiescent case to 50 (im for the 
lOOxBurst. For these wavelengths, the lOssenkopf &: Hennina (| 1994T ) dust opacities arc k ~ 100—400 cm 2 g" 1 , 



resulting in a required column density of dust for optical depth unity ~ 0.0025 — 0.01 gcm~ 2 . Assuming 
a gas to dust mass ratio of 100, and taking the fiducial envelope conditions above, this suggests that the 
peak emission is produced at an effective photosphere radius of i? p h ~ 50 AU. For comparison, the diamond 
symbols in Figure [2] denote the location of the emission-weighted r = 1 surface within each computed 
envelope. 



2.2. Timescales 

Three important timescales must be considered in order to determine the observability of protostellar 
variability from deeply embedded sources: the timescalc of the embedded source luminosity variations, 
expected to couple directly to changes in the accretion; the timescale for heating the envelope; and the 
differential photon travel times within the envelope. The timescale for variability in the luminosity of the 
underlying source is perhaps the most important feature. Large increases in luminosity can occur quickly in 
EXOr and FUOr outbursts. However, given our lack of a detailed model for such changes, we assume the 
simplest approach, that the central luminosity increases instantaneously to a burst level and then follow the 
effects of both the dust heating and photon travel times. Together, these processes work to spread out in 
time the observed luminosity jump, thus setting constraints on the types of luminosity variations capable of 
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being observed. We will return to the issue of the timescale for the embedded source luminosity variations in 
§[3] where we discuss the feasibility of monitoring the brightness of deeply embedded protostars as a means 
to discern the accretion history. 



2.2.1. Heating Times 

Figure [T] shows the required change in the dust temperature as a function of position within the pro- 
tostellar envelope when the luminosity of the central source increases, as expected during major accretion 
events. This temperature change, however, is not instantaneous as the dust retains a non-zero heat capacity 
and must absorb photons in order to heat up and re-emit in steady-state at the new temperature. 

The timescale for heating the envelope, as a function of radius, can be effectively broken into two parts. 
First, in the highly optically thick part of the envelope (i.e. when r < i? p h sec § I2.1[) a wave of energy will 
spread outward as fast as the dust can be heated to the new equilibrium temperature. Second, in the outer 
optically thin part of the envelope the heating time for the dust will be determined by a combination of the 
enhanced radiation field and the absorptive properties of the dust. 



Optically Thick Zone 



Ignoring for the moment the finite speed of light, the enhanced luminosity from the central source will 
propagate out from the centre at a rate fixed by the requirement for energy to be absorbed by the dust in 
order heat the dust to the new equilibrium temperature. The amount of energy deposited depends explicitly 
on the change in the temperature required, the amount of dust involved, and the heat capacity of the dust. 
Formally, 

£ abs (r)= / Airr 2 p(r)xC d (T)AT(r)dr, (8) 

where E a \, s (r) is the amount of energy absorbed by the dust out to radius r, Cd{T) is the specific heat capacity 
of the dust per unit mass, x = 0.01 is the fraction of envelope mass in dust, and AT is the required change 
in temperature. The lower bound on the integral, r denotes the distance from the protostar where the 
envelope begins (or alternatively the position within the envelope where the dust ceases to be evaporated). 
For the present analysis we do not vary rn in response to changes in the luminosity from the central source. 

For our fiducial m odel we have assumed that rp ~ 10 AU under the expectation that angular momentum 
or magnetic fields (e.g.. lTerebev et al ] ll984HGalli fc Shulll993h will produce a large degree of flattening in the 
inner region of the envelope. Observations of some young deeply embedded protostars sugg est that circum 



stellar disks may form early, although it is u nclear if these disks are rotationally supported (jJorgensen et al 
20091: IPersson et al.lEoii iTobin et al.ll2012l) . 



The dust heat capa city depends quite strongly on the specific properties of the dust assumed (jGuhathakurta fc Draine 



1989l : IChase et al.lll985h . however, a simple analytic form that fits most dust to within an order of magnitude 
for temperatures below ~ 500 K, sufficient as it will be shown for this analysis, is 



C d (T) = 10 



T 

Ik 



crg g 



K" 



(9) 



Given the strong dependence of the dust heat capacity on temperature, the heating time within the optically 
thick part of the envelope will be dominated by the most inner parts (i.e. r > r ). 
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To gauge the importance of the heating time, we now consider the amount of energy required to bring 
the inner few tens of AU in our fiducial model up to the equilibrium temperature, from zero, in the quiescent 
phase (see Figure[T]). The density at 10AU is pi(10AU) ~ 1.1 x 10~ 15 gcm~ 3 and the envelope mass out to 
10 AU is ~ 10~ 5 M© (this mass is contained in a fairly thin shell since r ~ ro), the total amount of dust is one 
percent of this value, and the dust temperature is ~ 100 K (see Figure[T]). Thus, E a b s (10 AU) ~ 2 x 10 35 erg, 
and the time to heat 

th = E a bs/L ps (10) 

is tfr,(10 AU) < 10 2 s. At 100 AU scales the envelope mass will be about 30 times larger, as we are still in 
the infalling portion of the envelope where M cnv oc r 3 ' 2 , but the dust temperature required for equilibrium 
will be about 3 times smaller (see Figure [T]) and the required heating per unit mass will be about 100 times 
smaller. Thus, as predicted, the heating of the optically thick region will be dominated by the inner region. 
For our fiducial model the time to heat the entire optically thick region will be th(R p h) 10 2 s. 

It is insightful to investigate the change in the envelope temperature and heating time produced by 
luminous bursts. The increase in the central luminosity manifests itself in an increase in the equilibrium 
dust temperature (see Figure [1]). However, the increase in temperature is almost exactly T oc L 1 / 4 such that 
the heating time = E a \, s /L ps remains effectively fixed. Figure [3] marks with circles the cumulative dust 
heating time for the optically thick portion of the computed envelope. 

We can also estimate the effect on the heating time produced by changing tq, the effective inner edge of 
the spherical infalling envelope. Given that T varies approximately as r -1 / 2 at small radii and p oc r~ 3 / 2 in 
the infalling zone, the integrand in Eqn. [8] is approximately proportional to r -3 / 2 . Thus i? a bs oc r 1//2 and 
varying the inner radius of the envelope by an order of magnitude, ro = 1 AU, only increases the heating 
time (Eqn. [TUJ) by a factor of 3. This, however, is likely to be an over-estimate of the increase in heating 
time with respect to ro since th e heat capacity used in Eqn. [9 i s a strong function of temperature only when 
the dust is below about 500 K ( Guhathakurta fc Drainel 1989 : Chase et al. 19851 ). By 1 AU the equilibrium 



dust temperature has reached this value even in the Quiescent case. Moving the inner edge of the spherical 
envelope in even further will therefore produce only a very modest increase in the heating time. 

Finally, it is useful to consider the time required to heat the gas as well as the dust, and the coupling 
time between the two components. The heat capacity of the gas is approximately 

C g ~ k/rrip, (11) 

where k is Boltzmann's constant and m p is the mean molecular mass. Thus, C g ~ 10 ergg K . At 
10 AU the equilibrium dust temperature is ~ 100 K and the gas heat capacity is about 10 times the dust 
heat capacity, per unit mass. Taking into account that the gas mass is 100 times the dust mass, the time 
to heat the gas in the envelope on this scale is about 1000 times longer than that required to heat the dust, 
t g ~ 10 5 s, assuming there is an effective way to couple the gas to the radiation field. 

The gas to dust coupling is most likely to take place via collisions between the dust and gas. The time 
for a dust particle to collide with molecular hydrogen is 

icon = {nadCsY 1 , (12) 

where n is the number density of molecular hydrogen, er^ = 7r r\ is the dust cross section, and r<2 is the 
typical dust radius. With each collision, the dust particle will lose energy to the molecular hydrogen. Given 
the factor of ten difference in heat capacities per unit mass, the dust particle will need to interact with an 
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equivalent mass of molecular hydrogen equal to ten percent of its own mass in order for it to appreciably 
cool. Therefore, the required number of collisions is 



where pd — 2.5 g cm 3 is the mean density of the dust. Together, the dust cooling time due to collisions with 
the gas is t c = t coU N coU : 

t c = 0.l( 4 ^) r A (14) 
\3piJ c s 

where pi = nm p . For pj(10AU) ~ 10~ 15 gcm~ 3 in our hducial model, and assuming r d ~ 10 _4 cm, the 
gas cooling time for the dust is t c (10AU) ~ 10 6 s. Given that this is much longer than the dust-photon 
coupling time, there is an effective separation of events, with the dust heating first and the gas reaching an 
equilibrium temperature much later. 

Putting all of the timescales together, the optically thick region of the envelope around the deeply 
embedded protostar should reach an equilibrium dust temperature in th(R p h) J$ 10 2 s, ignoring the finite 
speed of light which we return to in Section ^2.2.21 The gas will heat over a much longer time, t g (R P h) ~ 10 5 s, 
as it slowly draws energy from the already heated dust. 



Optically Thin Zone 

In the outer, optically thin part of the envelope the heating times are determined locally by comparing the 
energy input required to heat against the luminosity absorbed. Thus, for a shell at radius r in the envelope, 
with thickness dr, the amount of energy that must be absorbed in order to heat the dust is 

£abs(r) =4nr 2 p(r)xC d (T)AT(r)dr. (15) 

Equivalently, the energy absorbed per unit time from the interior by the shell is 

iabs(r) = Ltotdr, (16) 

where dr = pxndr and K is the absorptivity of the envelope per unit dust mass, averaged over all intensity- 
weighted wavelengths. The heating time th = -E a bs/-£abs is therefore 

th = ^C d (T)AT(r) t 

Ltot K 

which is independent of both the density p(r) and the shell thickness dr. 

Figure [D shows that most of the proto stellar luminosity is re-emitted near 50 pm at the effective pho- 
tosphere. The lOssenkopf fc Hennina (|1994l) opacity at these wavelengths is k ~ 400 cm 2 g -1 . Therefore, at 



a typical radial location in the envelope r ~ 10 3 AU, the dust temperature is ~ 10 K, and the heating time 
is t/j(10 3 AU) ~ 3 x 10 2 s. In Figure [3] we calculate directly this heating time from the calculated envelope 
properties as a function of position. 

The radial dependence for a centrally heated envelope shows th oc r 2 T A which is almost linear with 
radius since T varies approximately as r -1 / 4 at large distances. Thus, at 10 4 AU, the heating time for the 
fiducial model is ^(10 4 AU) ~3x 10 3 s (see also Figure [3]). Furthermore, during a burst the temperature 
in the envelope increases roughly as T oc L 1 / 4 . Therefore, as in the optically thick case above, the time for 
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the dust to reach equilibrium in the outer envelope is almost independent of the internal luminosity burst 
strength. 

In this outer envelope environment it is also useful to consider the time required to heat the gas as well 
as the dust, and the coupling time between the two components. For these low temperatures the difference 
in the heat capacity of the dust and the gas, per unit mass, is considerable, C g /Cd ~ 10 4 (see Eqns.|9]and 

At 10 3 AU, the collision time between the dust and gas, however, is much longer than in the extremely 
dense interior as the gas density is lower by a factor of ~ 10 3 . Compared with the interior, the time between 
collisions increases by ~ 10 3 due to the lower gas density while the number of collisions required to transfer 
significant heat from the dust to the gas decreases by ~ 10 3 due to the much lower heat capacity for dust at 
10 K. As such the dust cooling time due to collisions with the gas remains the same; £ c (10 3 AU) ~ 10 6 s. 

The time to heat the gas to the same temperature as the dust, however, increases dramatically. Com- 
pared with the dust heating time, the gas heating time is set by both the ratio of the heat capacity per unit 
mass and the ratio of total gas mass to dust mass: 

*» =arl (iO (18) 

Substituting in the fiducial values above, t g (10 3 AU) ~ 10 8 — 10 9 s. 

Putting all of this together, the optically thin outer region (r ~ 10 3 AU) of the envelope around the 
deeply embedded protostar should reach an equilibrium dust temperature on a timescale of 10 2 — 10 3 s, 
ignoring the finite speed of light (see H2.2.2[) . The gas will heat over a much longer time, ~ 10 8 — 10 9 s. 
Even in steady state, thermal balance calculations indicat e that the gas temper ature is decoupled from and 



lower than the dust temperature in outer regions of cores (jDoty fc Neufeldlll997l ). Nevertheless, the extreme 
separation of timescales for dust and gas heating seen here may have significant observational signatures (see 



2.2.2. Photon Travel Times 

In the above analysis of the heating time, we explicitly ignored the finite speed of light. The electro- 
magnetic propagation was assumed to be instantaneous, so that the heating time depended only on cither 
the time to indirectly heat the dust via contact with the expanding wave of energy (optically thick deep 
interior) or directly heat the the dust via absorption of photons emitted from the effective photosphere (outer 
optically thin envelope) . In this section we consider the importance of the finite speed of light in setting the 
timescales for observability of temperature variations in the envelope. 

When a change occurs in the luminosity of the deeply embedded central source, the increased photon 
flux travels outward through the envelope at the speed of light. The speed with which this wave traverses 
the deep, optically thick interior, will be set by the slower of either the speed with which the envelope heats 
or the speed of light. The dust heating time is ~ 10 2 s at 10 AU (see §2.2.1|) . implying a propagation speed 
of Wprop ~ 1.5 x 10 12 cms~ 1 , or fifty times faster than the speed of light. Thus, the actual progression of 
the photon flux, and the heating of the interior, should be set by the finite photon travel time. During the 
period when the interior optically thick part of the envelope is being heated, however, there is no observational 
evidence of the heating. Only when the heating wave reaches the effective photosphere are the reprocessed 
photons able to escape the envelope without being reabsorbed. 

Once the wave of heating reaches the effective photosphere, the radius at which the majority of the 
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photons produced locally arc capable of escaping through the outer envelope, the observed SED will begin to 
evolve in time. For an observer at a large distance from the source, this SED variation will begin at a fixed 
later time set by the arrival of photons from the near-side of the effective photosphere. For the observer, 
however, the entire effective photosphere brightness will not change instantaneously; rather, those parts of 
the photosphere further from the line of sight axis will take longer to arrive due to the additional geometric 
path which these photons must travel and the finite speed of light (see Figure 0] for a schematic diagram). 
For locations on this surface with an off axis impact parameter < i? p h, the time lag will be 



At(R b 



Rph 



1/2 



(19) 



The time over which the entire effective photosphere will be observed to heat up, and increase in 
luminosity, will be At ~ — ^ , or, assuming the photospheric radius is ~ 100 AU, At ~ 5 x 10 4 s. Thus, 
for an instantaneous increase in the central source luminosity the peak of the SED, which is dominated by 
emission near the effective photosphere, should gradually increase from pre-burst to burst conditions over a 
period of about half a day. 

Given that the dust heating times in the outer envelope remain fast, § I2.2.1[ the observed time response 
to the heating of the outer envelope will also be dominated by geometric effects. In general, the change in 
the emission of outer envelope shells should be spread in time by At ~ where the additional factor of 
two accounts for the fact that photons from the far hemisphere of the optically thin envelope are also able to 
reach the observer (see FigureH]). Since the outer shells at distances of 10 3 — 10 4 AU contribute significantly 
to longer (millimetre) wavelength emission, the timescale over which the SED evolves at these wavelengths 
will be weeks to months. 



3. Discussion and Complexities 
3.1. Observing a Burst 



As described in §2, an instantaneous increase in the central luminosity of a deeply embedded protostar 
will appear spread out in time to a distant observer primarily due to the geometric travel time of photons 
through the radially heating envelope (Figure EJ . The first evidence of the heating will be an increase in 
the luminosity of the source at wavelengths near, and shorter than, the peak of the SED (A ~ 50 /im). At 
these wavelengths, the SED will increase from its pre-burst equilibrium to a burst equilibrium in a matter of 
hours. At longer wavelengths the time for the SED to reach the burst equilibrium value will be significantly 
delayed, weeks to months, owing to the larger geometric paths for photons in the cooler outer envelope. 

Figure [5] reveals that the largest fractional change in the SED due to an increase in the central source 
luminosity occurs near the peak wavelength, where the strength of the variation is nearly identical to the 
variation in the burst luminosity. At longer wavelengths the difference between the pre-burst and burst 
SED is much less pronounced. In fact, for realistic, externally heated outer envelopes the evolution at long 
wavelengths is expected to be even less than that shown by the figure, since only internal heating of the 
envelope has been taken into account using DUSTY. 

The outer portions of real protostellar envelopes are heated by the interstellar radiation field and are 
found to have effecti ve temperatures betwe en 10 and 20 K ( Evans et al.ll2001 : Rosolowskv et al 1 l2008h . similar 
to starless cores (see lj0rgensen et alJl2006L for a discussion of why this is the case). As can be seen in Figure 
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[T] where again only internal heating is included, the majority of the heating of the very outer envelope even 
in burst scenarios is dominated by the external radiation field. This implies that the observable changes in 
the SED at long wavelengths will be smaller than suggested in Figured] 

Given the finite propagation of photons through the envelope, the observed SED as a function of time 
requires a proper time-dependent radiative transfer calculation in order to account for the variation in 
emission as the envelope heats. A much simpler approximation to the SED, however, is possible by assuming 
that the photons of interest are emitted primarily beyond the effective photosphere where the envelope is 
optically thin and that the timescale over which the heating occurs is short compared to the geometric travel 
times and can thus be assumed instantaneous (with a fixed delay). The SED can then be estimated simply 
by totalling the flux emitted by each radial and azimuthal envelope zone and by assuming the attenuation 
is negligible. This approximation does an excellent job at reproducing the longer wavelengths where the 
opacity throughout the envelope is extremely low, but overestimates the strength of the emission near the 
peak of the SED, since this emission is emitted near the effective photosphere and should be somewhat 
attenuated by dust in the outer envelope. 

Figures [5] and [5] highlight the temporal evolution of the SED for luminosity changes from Quiescent to 
10 x Burst and 100 x Burst, respectively, utilizing the simplified radiative transfer model described above. In 
each plot, the initial and final equilibria SEDs are presented, assuming no external heating of the envelope 
and thus correspond exactly to the SEDs shown in Figure [2] Over-plotted are observed SEDs as a function of 
time, where the zero-point is determined by the time at which photons from the burst first escaped through 
the effective photosphere and reached the Earth. For the plotted time-dependent SEDs, the outer envelope 
temperature is fixed at 10K in order to more accurately reflect real protostcllar envelopes. Thus, even at 
early times enhanced long wavelength emission is observed, compared with the D USTY models, because of 
the change imposed on the outer envelope temperature structure. 

Figures [5] and |6] reveal the quick increase in envelope luminosity at wavelengths near the peak of the 
see| and the lengthened time to reach the burst equilibria SED at longer wavelengths. The magnitude 
of the intensity change can be read directly from the y-axis labelling. In the case of the 10 x Burst, the 
fractional change in the intensity at mid-IR wavelengths (A ~ 50 fjxa) is ~ 10 and occurs on a timescale of 
a few hours. At far-IR wavelengths (A ~ 200 (im) the fractional change in intensity is ~ 5 and occurs on 
a timescale of ten days. Finally, at sub-mm wavelengths (A ~ 500 fxxxx) the fractional change in intensity is 
~ 2 and occurs on a timescale of a few months. For the 100 x Burst case, the fractional intensity changes 
increase significantly, while the timescales over which the changes are observed remain roughly constant 
(as anticipated in the previous section). We note that near- and mid-IR wavelength observations might be 
challenging because m uch of the detected emission is produced in outflows rather than in the envelope (e.g. 
J0rgensen et al. 20091 ) . which could confuse any outburst signature. 



The time-dependent SEDs shown in Figures [5] and [5] require that the luminosity burst of the deeply 
embedded protostar persist over the observed time period. If the burst is short in duration then a wave of 
heating followed by a wave of cooling will propagate through the envelope and the observed temperature 
in any radial and azimuthal envelope zone will be much more complex. Such time-dependent models are 
extremely interesting but beyond the scope of this paper. 



2 The simplification to the SED described in the previous paragraph is responsible for the strange behaviour observed at 
the shortest wavelengths. In particular, in the simplified radiative transfer model the shortest wavelengths reach intensities 
somewhat larger than expected for the burst equilibria. This occurs precisely because the attenuation of these photons through 
the outer envelope is not taken into account. 
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3.2. Comparing Dust and Gas Temperatures in the Envelope 



As discussed in detail in t|2.2.1[ the dust responds much quicker than the gas does to a burst in luminosity. 
For the inner optically thick zone of the fiducial envelope model, the lag in heating the gas (At ~ 10 5 s) is 
offset by the photon travel time (At ~5x 10 4 s) and observing the decoupling of the gas and dust is likely to 
be difficult. For the outer envelope, the optically thin zone beyond the effective photosphere, the gas heating 
time can take many years, significantly longer than the photon travel time (months). Thus, depending on 
the cadence of the underlying luminosity variations, and despite the relatively high gas density p s > 10 5 cm 3 , 
it is entirely possible that the outer envelope dust and gas may not attain temperature equilibrium. 



Kim et al.1 (j 2 1 lh and IVisser fc Berginl (|2012al ) attempt to characterize the signatures of past outbursts 
in the chemistry of envelopes. While some excitation o f the outer disk ma y lag the outburst by years, 
many lines are produced by outflow-envelope interaction (jVisser et al.ll2012bl ). where higher densities may 
improve the coupling between dust and gas and photochemistry would occur quickly because of the presence 
of outflow cavities. 

Although beyond the scope of this paper, it will be extremely illuminating to connect the timescales 
and strengths of protostellar l uminosity variations due to bursts as suggested by specific theoretical models 
for protostellar evolution (e.g. lDunham fc Vorobvo vl l2012h with the time dependent changes in both the dust 
and gas temperatures in the envelope. Such investigations could be used to define survey properties such as 
target cadence, continuum wavelength coverage, and which molecular lines to observe. Alternatively, a dust 
continuum and molecular line monitoring campaign of known deeply embedded protostars would inform the 
theoretical models (see <j3. 41 below). 



3.3. Alternative Envelope Models 

The fiducial model presented above is extremely simplistic, relying on spherical symmetry and power- 
law radial density profiles. The results, however, are much more robust than might be initially assumed. 
First, retaining the inside-out collapse model and changing only the effective temperature produces little 
change in the time dependent solution for the SED since the timescale is set primarily by the photon delay 
as a function of radius in the envelope (see Eqn. [T9]) . 

For the inner, optically thick envelope this situation will change only if the mass and equilibrium 
temperature beneath the effective photosphere becomes extreme enough to require a significant increase in 
the heating time (see Eqns.l51and [TU|) . Given the expected range of temperatures defining low-mass protostar 
envelopes, T ~ 10 — 20 K, this is unlikely to occur (see Eqns. Q] and [2]). In the outer, optically thin region 
of the envelope, the effect of changing the envelope structure is diminished even further because the only 
relevant terms are the envelope temperature and the protostellar luminosity (see Eqn. I17[) , and those terms 
tend to counteract one another. 

The situation for deeply embedded massive protostars is less clear, precisely because the physical con- 
ditions of these objects are not yet well understood. If, however, the effective temperatures are significantly 
higher than for low-mass analogues then the density in the inner regions will be much larger and the heating 
time could, in principle, become longer than the photon delay time within the optically thick zone. As well, 
the effective photosphere is likely to increase in size due to the enhanced envelope opacity (directly related to 
the density) and this will increase the time over which the peak of the SED brightens, Nevertheless, as noted 
above, the heating time in the outer envelope is less likely to be strongly affected by either the increase in 
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envelope density or equilibrium temperature and the results above, perhaps scaled to larger envelope radii, 
should still hold. 

Considering again low- mass protostar s, the above commentary holds for more realistic envelope profiles, 
such as collapsing Bonnor-Ebert spheres ( Bonnor 19561 : Ebert 1955 : Myers 2005 ). since the enclosed mass 
and equilibrium dust temperatures are not significantly different than for the simpler singular isothermal 
sphere model. Nevertheless, while beyond the scope of this present paper, it would be useful to produce 
detailed time-dependent SED profiles for these often used envelope models in order to measure the finer 
distinctions and deviations. 

Edge on disks around these deeply embedded sources will further enshroud the central protostar but will 
not be particularly important in determining the heating of the envelope or the observability of the changing 
SED unless the disk is extremely large, -R disk ~ Rph- Alternatively, outflows driven by accretion p rocesses 



carve out cavities inside the envelope (e.g. lSnell et al.lll984t |j0rgensen et alj|2007t iTobin et al.ll2010l ). When 
the outflow axis is oriented with our line of sight, direct evidence for protostellar variability may be detected 
quickly. For envelopes with outflow cavities unaligned with the observer, scattering of short wavelength IR 
photons off the inner wall of the outflow cavity may also provide an early beacon of burst activity. Since 
the only delay in observing the beginning of a burst through an optically thick envelope is due to photon 
propagation, the major advantage of seeing directly to the protostar is that the observed change in intensity 
at short IR wavelengths should be almost instantaneous, dependent only on the size of th e emitting regio n 



around the protostar and expected to be ~ R ps . Some FUOr objects still retain envelopes (jZhu et al.ll2008f >. 
which must be optically-thin in our line of sight - perhaps because of an outflow cavity. 

Finally, the time-dependent SED could in principle be inverted to measure directly the density structure 
within an observed envelope, independent of theoretical models. In theory, combining time-sampled SEDs 
across a broad range of wavelengths should allow for a measure of the emission variation associated with the 
parabolic region trace d out by the time-delay (see Figure IH). Reverbe ration mapping has been successfully 
applied to AGN (e.g. Blandford fc McKeei 1982 : Peterson et al. 2004 ) but has not yet been applied to en- 
velopes or disks around young stars. Given that the expected equilibrium temperature profile in the envelope 
depe nds predominantly on the d ust opacity and the embedded protostar luminosity, assuming the former 
(e.g. lOssenkopf fc Henninglll994l ) and measuring the latter should provide leverage on such an inversion. It 
remains uncertain as to whether the required sensitivity to these differential measurements is within the 
limitations of current instruments and below the non-equilibrium variations expected within the protostellar 
sources themselves. 



3.4. Survey Possibilities 

The results present in Sj2] and discussed in more detail in §3.1| 13.21 and 13.31 suggest that a long-term 
monitoring campaign of deeply embedded protostars would constrain the possible luminosity variations, and 
hence mass accretion rates. Indeed, while the low luminosities observed for most deeply embedded protostars 
suggests that mass accretion is episodic, arriving as short bursts with extreme accretion rates, very little is 
actually known about the cadence or the relative number of bursts with varying mass accretion amplitude. 

Thus, while it is possible that deeply embedded protostars are quiescent and non-varying for most of 
their evolution, with very occasional strong bursts in which the mass is deposited onto the protostar, it is 
equally likely that these sources are constantly varying at low mass accretion amplitudes. With present 
instrumentation it should be relatively easy to monitor these sources for fractional intensity variations at or 
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abovc fifty percent (and likely as low as ten percent), providing significant leverage on low amplitude mass 
accretion variation. Here we suggest four instruments that might provide observations which would help 
discriminate between the above possibilities. 

Space-based imaging with the Herschel Space Observatory (Hcrschel), now in its last months of obser- 
vations, has produced an impressiv e archive of star-forming molecular clo ud observations using both the 



PACS and SPIRE instruments (e.g.. lAndre et alj|2010t iMolinari et al.ll2010h - imaging at wavelengths near 



the expected peak in the observed SEDs of deeply embedded protostars, where the expected variability 
has the largest intensity amplitude, as well as at longer wavelengths. In many cases multiple observations 



of th ese protostars have been undertaken and already variability at 70^m has been observed ([Billot et al 



20121 ). While it is unfortunate that Herschel will not be available for a carefully defined survey in search of 
protostellar variability, the archive should contain a trove of useful observations for future data-mining, In 
particular, sources shown to significantly vary across Herschel observations should be specifically targeted 
for follow-up with ground-based observations. 

Observing variability may be possible from the ground, but only in the driest locations. The James 
Clerk Maxwell Telescope with the SCUBA-2 camera is capable of surveying lar ge (~ 43 square arc min) 
regions within molecular clouds at 450 fim and 850 /jm with excellent sensiti vity ( Holland et al. 20061). A s 



part of the telescope legacy, a large Gould Belt Survey is already underway (jWard-Thomspon et al.l 120071) . 
however, there are as yet no plans for systematic time-sampled observations of deeply embedded protostars. 
A dedicated survey of a few small clusters of deeply embedded protostars should be both efficient and effective 
at uncovering variability, as each observation should take only a few minutes and differential photometry 
should simplify the data reduction. Similarly, the 25m CCAT telescope under development is planned as a 
single-dish telescope at 5600 m on Cerro Chajnantor with first light instruments expected to include wide- 
field sub-mm cameras operating at < 350 /im. Since the relevant wavelengths are longward of the peak SED 
of envelopes, monitoring with SCUBA-2 or CCAT should be separated in time by weeks (see Section l3~T]) . In 
the event that CCAT's camera allows 200 /im observations, the cadence of monitoring should become almost 
daily. 

Finally, the ALMA Observatory provides interfcrometric sub-mm observations and is capable of imaging 
the effective photosphere of nearby deeply embedded protostars. At 150 pc distance, the effective photo- 
spheres of the nearest embedded protostars will be ~ 1", well within the range of spatial scales captured by 
ALMA. ALMA observations will therefore enable wavelength dependent, spatial dependent, and time de- 
pendent observations of protostellar variability. Additionally, ALMA can monitor both the dust continuum 
and individual molecular lines simultaneously, allowing for an investigation into possible deviations between 
the dust and gas temperature. While using ALMA as a survey instrument to find accretion variability is 
unlikely, the power of ALMA could be utilized extremely effectively in following up bursts found by the 
JCMT or CCAT. 



4. Conclusions 

The growth of protostars may occur stochastically in large outbursts, while even in quiescent states the 
weak accretion should vary with time. In this paper, we investigate whether and how to detect changes in 
accretion onto protostars. When the central source increases in luminosity, the envelope brightens at all 
wavelengths and the SED peak shifts to shorter wavelengths. The dust quickly reaches this equilibrium in 
both the inner and outer envelopes. The timescales for changes in the SED are thus determined by the photon 
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travel time and how quickly energy can propagate outward. These changes would be easiest to detect in the 
mid- to far-IR, occurring on timescales of hours to days, but should also be possible at sub-mm wavelengths, 
occurring on timescales of days to months, with ALMA, the JCMT or CCAT. The gas, however, reaches 
equilibrium much more slowly due to the slowness of the gas to dust coupling and the much larger heat 
capacity, per unit mass, of the gas. 
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Tabic 1: Physical properties of the fiducial envelope model 



Property Value Reference 

T (K) 10 initial condition 

t x (yr) 1.6 x 10 5 initial condition 

M onv (M Q ) 1.5 initial condition 

M ps (M Q ) 0.25 EqnH 

L ps (L Q ) 1.2 Stahler (1983) 

L acc (£ ) 5 EqnE] 

ilOXBurst [Lq) 12 

-^lOOXBurst (£©) 120 

R cnv (AU) 1.8 x 10 4 Eqn[5] 

(AU) 6 x 10 3 Eqn[3] 

R ph (AU) -50 t ~ 1 



^ (r) (gem" 3 ) 1.1 x 10- 15 (r/10AU) 3/2 Eqn[U 
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Fig. 1. — Dust temperature as a function of radius for our fiducial envelope (see Table [TJ as given by the 
radiative transfer code DUSTY (see text). The quiescent case (blue solid line) assumes L ps = 1.2L and no 
accretion luminosity. The two different burst scenarios, lOxBurst (red dashed line) and lOOxBurst (red dotted 
line), assume Ltot = 12 and 120 L©, respectively. The two solid black vertical lines indicate the (externally 
fixed) inner evaporation radius of dust and the outer edge of the envelope, whereas the dashed vertical line 
marks the boundary between the static and infalling density regimes (see Eqn[3]). The effective photosphere 
of the envelope is shown as a diamond marker for each model and moves outwards with increasing central 
luminosity. Note that DUSTY calculates the temperature solely from internal heating, without an external 
radiation field. The red horizontal line indicates the lower limit to the envelope temperature which we will 
later apply to calculate time-dependent SEDs (Figures [S] and IB]) • 
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Fig. 2. — Equilibrium Spectral Energy Distributions as obtained from DUSTY for the envelope properties 
shown in Figure [TJ Around the peak wavelength the brightness varies almost linearly with the total source 
luminosity, Ltot- At long wavelengths the brightness varies significantly less (see text). 
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Fig. 3. — Heating time for the envelope as a function of radius. The red open (solid) dot indicates the 
integrated heating time out to the effective photosphere, denoted by the left vertical dashed line, for the 
lOxBurst (lOOxBurst) scenario. The red dotted (dashed) line shows the individual heating time of the 
optically thin outer shells as a function of radius for the lOxBurst (lOOxBurst scenario). The drop in 
heating time at the outermost radii is due to the imposition of the minimum dust temperature of 10 K which 
significantly reduces the temperature differential between the pre- and post-burst at the largest distances 
from the protostar (see Figure [T]). 
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Fig. 4. — Schematic visualization of the observed heating process. Due to significant differential photon 
travel times as a function of radius and inclination angle in the envelope, the observed heating of the envelope 
is geometrically distorted (see Eqn [T9| . The grey shading indicates times of Oh (black), lOh, lOOh, lOOOh, 
3000h, and 6000h (white). Time zero is defined when the first photon reaches the observer, located in a 
positive x direction in the above schematic. The time delay of photons is the dominant factor in determining 
the time-dependent SED following a burst (see Figures [5] and HJ) . 
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Fig. 5. — Time variation in the observed SED of a deeply embedded protostar after a lOxBurst increase 
in the source luminosity (see text). The dotted lines show the respective equilibrium SED for the quiescent 
case (solid blue) and the lOxBurst (empty red) as obtained with DUSTY. The grey shading follows the 
change in the SED over time as calculated assuming an idealized envelope, optically thin at all wavelength 
(see text). This assumption produces almost identical results with DUSTY for wavelengths longer than 
~ 100 /im (with the exception of long wavelengths in the Quiescent state - see below). For significantly 
shorter wavelengths the optically thin assumption is incorrect and the amount of emission is overestimated. 
The shading corresponds to to 0.5h (dark grey), 5h, 50h, 500h, and 5000h. After fifty hours the flux around 
the peak wavelength at near 100 [im has already changed by a factor of 10, whereas there has been little 
change at significantly longer wavelengths. The long wavelength region of the spectrum which experiences 
an extra lift, compared to the DUSTY SEDs, due to the assumption of external heating and a minimum 
dust temperature of 10 K is hatched for clarity. 
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Fig. 6. — Same as Figure [5] but for the case of a lOOxBurst. 



